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The disturbances produced in a homogeneous plasma by passing charges is analyzed. 

1. Introduction 

A source which is effective during a limited period can be described by a function of time 
or by an equivalent Fourier spectrum. It is well known that the disturbances produced by 
such a source depend mainly on its highest frequencies shortly after the arrival (determined 
by a propagation velocity) of the first effect, at t=t say, whereas the lower frequencies de- 
termine the effects observed later on. The Fourier synthesis of the disturbance, which has to 
include the contributions of all frequencies, will represent the special effect of the low frequencies 
in an obscured form even when the latter finally dominate. A new representation is then 
desired; it may be obtained as follows. 

The time dependence of the disturbance may be known as a Laplace integral 
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which is equivalent to a Fourier integral with a path of integration L parallel to the real axis 
of the a>-plane. This integral can then be reduced in two different ways: 

(a) a parallel shift of L to the real axis itself. This results in a conventional Fourier 
synthesis in terms of real frequencies. In the case of poles on the real axis, these are to be 
avoided by an indentation just above it. In the limit of an infinitesimal indentation this 
reduction leads to the sum of a main contribution represented by a Oauchy's principal value, 
and of additional contributions consisting of half the residues at the poles. 

(b) a transformation of L into a closed contour which results by adding to L the upper 
half of the infinite circle of the co-plane if iK^o and its lower half if ^>/ . In the first case the 
integral will vanish, as it should, if /(co) has no singularities whatever above L. However, 
the other contour to be applied for t^>t , encloses the singularities of /(co) situated in the half 
plane Im co<0 or possibly on the real axis. The integral itself can then be reduced to the 
residues at the enclosed poles and to contributions along branch lines connected with the 
other singularities. If the branch lines do not extend up to infinity, we get by contraction of 
the integration path around the poles and branch lines under consideration contributions 
connected with discrete frequencies and with frequency bands respectively. 

As an example we consider the Laplace integral 
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A shift of the integration path to the real axis here results in: 
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All frequencies are contained in this integral, though the particular frequency co has a finite 
amplitude even when observed in an infinitesimal band around o> . On the other hand, trans- 
formation to a closed contour as suggested in method (b) reveals the vanishing of the integral 
for £<0, while it simply reduces to (—2tt/o)o) sin (a> £) for £>0. Thus the property of the 
signal of containing only a single frequency co right after its start at £=0 was hidden in the 
conventional Fourier representation (2). 

Idealized ionospheric disturbances, analyzed with the aid of the above procedure, will 
be discussed in the next sections. They refer to a single infinitesimal current element, only 
effective at a special moment, and also to a charge moving with constant speed along a straight 
line. The latter situation simulates a charged particle emitted by the sun, that travels along 
a magnetic line of force through the ionosphere. The limiting cases of a zero and of an infinite 
magnetic field are worked out in detail, the intermediate case of a finite magnetic field is only 
mentioned briefly. The ionospheric plasma is assumed throughout as homogeneous (constant 
plasma frequency ai pl ) and cold. 

2. Disturbance Generated by a Current Element in a Plasma Without Magnetic 

Field 

2.1. Derivation of Laplace Integrals for the Field 

Let the current element be situated at the origin, and directed along the 2-axis, while 
effective only at £=0. This corresponds, in the case of a normalized momentum 8(t), to a cur- 
rent-density distribution: 

I=8(x)5(y)8(z)8(t)Z 2 , (3) 

-> 

u z being the unit vector in the 2-direction. 

The effective diefectric constant e=l — <x> 2 pl /u> 2 for a homogeneous pfasma at frequency w 
then involves the following Maxwell equations in gaussian units: 

Curl E+ c ~M = °> 

,dH 1 / d 2 , , \ £ 4tt d/ 

=^8(x)5(y)5(z)8'(t)u,. (4) 

The complete solution in terms of Laplace integrals reads: 
E=—4ir8(3c)8(y)8(z) cos (w p it)U(t)v z 
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The sign of the square root is defined by the property -yjoo 2 — co^->co along the infinity circle, 
while a crosscut is assumed along the section — co^<co<w^ of the real axis. Further, u r and 
Up mark unit vectors in the radial and azimuthal directions corresponding to polar coordinates 
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r, 0, <p with z=r cos 0, x J riy=pe i<p =r sin 6e itp ) 

U(t) represents Heaviside's unit function, while 8 may be any positive number. 

The solution (5) can be verified as follows. We evaluate curl E (2) with the aid of the vector 
identity: 

curl curl curl = — curl V 2 , 
as well as the relation: 

ikr 



(V 2 +F)^- =4«r6(x)8(y)8(x). 



We thus obtain: 

curl E™=^- 2 curl/^ T^ 4r-W^- 2 -v\ 

u g 8(x)8(y)8(z) du — 2 r V- 
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The last integral equals —2iriU(i) cos (cjpit) from which the second term proves to reduce to 
— curl E {1 \ Hence we find for the total E field: 

curl E=-^- 2 curM -^ efoe-'-'w* V- 1 --* 1 >■ 

The first Maxwell equation of (4) is then readily verified by applying the further identity 

-> -> ► -» 

sin • X ' (r) • i^ = u z X grad X (r) = — curl { u z • X (r) } (7) 

to the magnetic field i?. 

The other Maxwell equation can be checked, without working out the curl curl operator 
-» 
occurring in curl dll/dt when (5b) is replaced by a curl equivalent to it according to (7) . In 

-> 
the further elementary derivation the application of the operator d 2 /dt 2 -\-o) 2 P i to E {2) amounts 

to a cancellation of the denominator of the integrand. 

The first contribution E {1) of the electric field is recognized as a nonpropagating field with 

a frequency equal to the plasma frequency; it can not move away from the source. The main 

-> -» 

contribution E (2) , as well as H, is composed of spherical waves of complex frequencies. 

2.2. Frequency Spectrum After the Start oi the Field 

An application of method (a) of the first section shows how all real frequencies contribute 
to the final field, but the explicit form of this spectrum is not very instructive. The method 
(b), however, directly shows the predominant role of the lower frequencies after the signal has 
started at the point of observation. In fact, the approximation 



for the exponential factor in the integrand of both integrals of (5) involves that the integration 

V T 

path can be closed along the upper or lower half of the infinite circle when t<Cr or £>- respec- 

c c 

tively. In the first case the only singularities of the integrands, situated at u=a) P i and co= — co pl , 
are outside the contour and the integrals vanish. In the second case the contour can be con- 
tracted to a loop L around the part — co^^w^co^ of the real axis which was introduced as the 
branch line fixing the sign of -yjco 2 —^; the ends of this line constitute the only other singulari- 

ties of the integrand, viz, its two poles at u v i and — a> P i. Hence the main contribution E i2 \ 
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e.g., can be represented at all times by: 



# (2) =?r-curl curl i 
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The contraction of L involves the sum of the complete residues at the poles, and of the 
limit for e-^+0 of an integral along the interval — co /jZ + e<co<co 7?z — e of the real axis; along 
the latter we have to take into account the new values i^u 2 pl — u 2 and — i^oo 2 pl — w 2 of -yJuP—upi 
just above and below the branch line. Both residues can be added together to a term — 2iri 
cos (a pit). The remaining integral can be reduced to an integration over the positive interval 
0<^oo<^to P i— e only, which, however, involves an integral yet to converge for e=0. The residues 

are missing in the corresponding integral for H. We thus arrive at the following final expres- 
sions without an e limit: 
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For £^>- we infer a dependence on the lower frequencies co<co pZ only. 

2.3. The Field Expressed in Terms of Bessel Functions 

The integrals in (5) are connected with Bessel functions. This can be shown with the aid 
of the operational relation [Van der Pol and Bremmer, 1955] 

JoUl^?)-U(t-r)=-^= 6-VFFI (Re g >0), (9) 

V2+1 

in which we use the nomenclature h(t)=f(q) for 

/(fl)=2 P dte-*h(t). 

We replace t by co p ^ [which requires the replacement of q by q/u P i] and also the nonoperationai 
variable r by co^r. We then obtain the following extension of (q): 
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(Rej>0). 



(10) 



We next apply to both sides the operator (likewise nonoperationai) dr, which leads to: 

J Tic 

u ( tJ i)L dTj ^^ II?) - q m^- (Re 2>0) - (11) 

The inversion integral of this new operational relation reads: 

2ttJ_ 00 + z5 (wji— co 2 ) 
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(12) 
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The derivative of this relation with respect to t yields at once a new expression for the E' 2) 
contribution of (5). We find: 

U (t--) r , 

' J r/c 



£ ,(2) =^ 7 curl curl 



y • 



(13a) 



A corresponding reduction for the magnetic field represented by (5b) is arrived at by first 
deriving from (12): 



(£+«&,) {u (*-=)£ dr • j o( ^,v?=?5 } =i £ + ; 



da • e 
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We next apply the operator — d/dr to (10) while taking t=-> and also the operator (d 2 /dt 2 

+ <4;)=(g 2 +o4 z ) to (11). In both cases we get the same operational "image" in the right- 
hand side which involves an identity for Bessel functions according to which (14) can be replaced 

Hence the magnetic field (5b) can be expressed as follows in terms of Bessel functions, without 
any integral: 



^^•■IIJI^W^^O-O}]-^ 



(13b) 



3. Disturbance Generated by a Charge Moving Along the z-Axis Through a 

Plasma Without Magnetic Field 

3.1. Derivation of Laplace Integrals for the Field 

We assume a unit charge moving with a constant velocity v along the 2-axis according to 
the equation z=vt. The corresponding current-density distribution is given by: 



~-d(x)d(y)d (*-£) v> t . 



The connection with the nonmoving current element of the previous section is shown by 
the equivalent representation: 






dt -6(x)5(y)8(z-vto)5(t-to)-u t 



Hence this new current distribution is obtained from the former (3) by first substituting z—vt 

Joo 
dt . In view of the linear dependence 

-> -> 

of the field vectors E and H on the sources, the fields for the present situation may be derived 

by applying the mentioned substitutions and operator to all results arrived at in section 2. 

We first consider the electric field. In the case of (5a) the mentioned procedure results 

in the following expression if the t integration is evaluated in the first contribution EQ), and 

if the order of the integrations with respect to co and t is inverted in the main contribution E( 2 ) : 



E=-47r8(x)d(y) cos ^ a> pl (*-") V ■ U (*—) " ^ 

r 
+h- curl curl ^ v z da> - 2 dt 



x« + ?-V^ 2 -co^2Vp24-(2-^o) 2 
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The substitution 



t = — i - sin < X+i arc sinh 
v v I 



/ 2 2~~ 

•yco — co p ; 



(5-0 
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reduces the ^-integral to a Sommerfeld integral for a Hankel function H$ L) (iz)=— K (z) which 
is given'explicitly by 



;^{>V(B)«'+*} 



provided 'that^the modulus of the argument of the square root is smaller than tt/2. This in- 
volves the*following final expression for the electric field: 



E=-4Td(x)8(y) cos < a, pl (t-f)\ ' U (t—f)' ^ 



H — curl curl 



U z 
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(16a) 



The corresponding reduction of the magnetic field can be carried through by first trans- 
forming the formula (5a) for this field in the case of a nonmoving current element into 



H=- — curl 

2irc 



{|j^rfa, e ~^V^^ 



which is equivalent to the original expression in view of (7) . The resulting expression for the 
magnetic field of the moving charge then becomes: 



in the derivation of which we have applied a vector identity similar to (7), viz: 



(16b) 



/(s)-0 7 (p) -^ = -curl {u z -j(z) -flf(p)}. 

3.2 Discussion of the Frequency Spectrum 

The integrals in (16) cannot be transformed into contour integrals, in contrast with the 
corresponding Laplace integrals in (5) for a nonmoving current element. Mathematically 
this follows from the exponential behavior at infinity of the part 



^•M^B)"^} 



of the integrands. In view of the asymptotic approximation of K the total exponential factor 
of the integrands can then be represented for |co|-h>oo by 



-* w (*-f)- p y^-2. M 



For v<^c this involves an infinite increase of the integral along those parts of the infinite circle 
of the co-plane at which 



tan (arg w )> p yI-I^--£). 
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Both the upper and lower half of the infinite circle therefore contain a sector in which the 
integrand increases exponentially. The physical reason for the impossibility of obtaining 
contour integrals which are different for large and small t is the lack of a starting time before 
which the signal vanishes. The situation would be quite different in the opposite case v^>c 
where the disturbances propagating with velocity c move behind the source; this latter situa- 
tion is realized for Cerenkov radiation. 

We shall further restrict ourselves to the case v^c. We then can only apply method (a) 
of the procedures mentioned in the introduction, which amounts to a shift of the integration 
path of (16) to the real axis, or to determining the limit for 8-^0. In (16a) we then have to 
take, in view of the poles at u=±u ph the sum of Cauchy's principal value (V.P.) connected 
with the latter, and of half the associated residues. Moreover, the path of integration can 
finally be reduced to the positive half of the real axis. In (16b) the poles are missing, and it 
is unnecessary to consider a principal value. The final results become: 
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Obviously, the first term of E represents a contribution left behind along the trajectory 
which has been covered by the moving charge up to its position at the moment of observation. 
All contributions only depend on p and t—zjv, which indicates that the field moves along with 
the same velocity as the charge. The spectral density of the finite contributions depends first 
of all on the monotonically decreasing K function. This decrease determines the predominance 
of the lower frequencies. As a matter of fact, only those frequencies contribute considerably 
for which the argument of K does not exceed the order of magnitude of unity. Working out 
the corresponding condition in the case v << c we infer that the field is mainly restricted to 
a cylinder of radius c/co pi around the trajectory of the charge (2-axis) while the cutoff fre- 
quency decreases rapidly inside this cylinder when moving away from its axis. Only very low 
frequencies may yet penetrate a little outside the mentioned cylinder, up to distances roughly 
of the order of v/o. 

3.3. The Field Expressed in Terms of Elementary Functions 

In the case of a nonmoving current element the field could be reduced to Bessel functions 
(see section 2.3) . For a moving charge the magnetic field can even be connected with a simple 
algebraic function by substituting in (16b), for 8— >0, the expression 



*W(B)^H/>-- ? '- w "* 



sinh <t> 



This relation is again connected with the representation of Hankel functions by a Sommerfeld 
integral. 

After the substitution we invert the order of integration with respect to co and <p. The 
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integral over co then amounts to a delta function according to which we find: 

f°° , - p -<* P i'Cosh4> / z /l 1 . , 



.!' 



zco • e 



This integral can be worked out most conveniently by introducing the argument of the delta 
function as a new integration variable, and applying next the sifting property of integrals con- 
taining such a function. The formula (16b) thus finally proves to be equivalent to 
the representation: 



S=V(i-i)-("0' 



-> 
H= 



c dp 



V(B>«+(<-0' 



The exponential factor, which is missing in the absence of a plasma, indicates that the 
latter compresses the field in space and in time. 

4. Disturbance Generated by a Current Element Along a Line of Force of 
an Infinite Magnetic Field Embedded in a Plasma 

4.1. Derivation of the Laplace Integrals for the Field 

According to the Appleton-Hartree formula the external magnetic field only enters into 
the refractive index of the plasma in the ratio go/go h (u h = electron gyrofrequency) . Hence, the 
limiting case of an infinite magnetic field is related to that for very low frequencies and therefore 
is instructive for these frequencies. 

Let the homogeneous external field be in the z direction, and the current element of nor- 
malized momentum at the origin. Its current-density distribution is then given, once again, 
by (3). The dielectric constant of the medium under consideration is given by the tensor: 
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This tensor involves the Maxwell equations: 
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curl E-\ ^— =0, 

c dt 
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The solution corresponding to (5) in the case without magnetic field reads: 
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(18a) 
(18b) 



(19) 
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The first equation (19) is verified readily, when representing the H field as a curl vector, 
applying (7). In order to check the second equation of (19), the curl curl operator then occur- 
ring in its first term may be split into curl curl= — V 2 +grad div. The Laplace operator 
V 2 operating on the quantity 



V.V2- 



VwV — co 2 pl p 2 
should next be transformed with the aid of the relation: 



d 2 1 d 2 



W + ^ 2+ 



<*?>z J Jx 



k Vz 2 +2/ 2 +«:2 



for 



-y +cxz 



■\a 



5(x)5(y)8(z) 



(argV«= 2 argaj 



(20) 



a=oj 2 /(o) 2 — co 2 pl ) and k=^/co 2 —co 2 pl /c. 



The term with delta functions then leads to the right-hand member of the second equation 
(18) whereas the other terms prove to cancel. Thus the correctness of (19) depends first of 
all on the validity of (20), a relation which can be proved, e.g., by integrating it (applying 
the method leading to Gauss' integral theorem) over a volume element around the origin; 
this relation constitutes an extension, for a^l, of the property (6) for the point-source solution 
of the three-dimensional Helmholtz equation. 

4.2. Frequency Spectrum After the Start of the Field 

The derivation is completely similar to that in the case of the absence of a magnetic field , 
used when deriving (S) from (5). The situation is even simpler since the branch point of 

(co; ;/ p 2 — or/- 2 ) 172 at co=co pl p/r 

constitutes the only singularity, as the pole is now missing even for the electric field. The 
results can be represented by: 

r /i o 
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^oo 2 plP 2 -co^~ 



(21) 



Just as in the case without magnetic field, described by (8), the field starts at the moment 
t=r/c, but the cutoff frequency now equals co pl sin 6 instead of co pi . The amplitude is maximal 
near this cutoff. This implies a concentration of the lowest frequencies to directions very 
near to that of the magnetic line of force passing through the radiating source. Here we infer 
a similarity to "whistler" propagation. 

4.3. The Field Expressed in Terms of Bessel Functions 

The expressions corresponding to (13) in the absence of the magnetic field are obtained 
here with the aid of the relation: 



da* 

-oo + iS 



V 'co 2 — co\ 
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of which (15) constitutes the derivative with respect to v. Applying this relation for co pZ 
replaced by u P i sin 0, we obtain from (21): 



M^s-tS) H-D 



Jo(«piSin0 ■\ t2 — 7 -l) 



>> 
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5. Disturbance Generated by a Charge Moving Along a Line of Force of an 
Infinite Magnetic Field Embedded in a Plasma 

5.1. Discussion of the Frequency Spectrum 

Here, again, the field for a moving charge can be deduced from that of the nonmoving cur- 
rent element, just as discussed in section 3.1 in the absence of a magnetic field. Starting from 
(19) we now have to evaluate the following integral over t \ 



L- 



■/: 



dt 



iwfo-h- V(w 2 -w P i 2 )p 2 +a>2(z-^o) 2 
V(to 2 -co^)P 2 +co 2 (3-^o) 2 



Again, this integral could be reduced with the aid of Sommerfeld integrals to Bessel and 
Hankel functions. However, it is impossible to get a single representation for all co-values of 
the integration path of (19) which, therefore, has to be split into different sections. This 
splitting is simplest in the limit 5^0 and can then be reduced to the positive part of the real co-axis. 
The final results read as follows in terms of Neumann and if functions with positive arguments: 
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-i £ <Z<o cos {co («-*)}* (' VH^^)}} (22) 
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These expressions can be verified most conveniently by showing that they do satisfy the 
Maxwell equations (18), with the right-hand side of (18b) replaced by 



^ 1 z b(x)5(y)&' (*-£)• 



(23) 



The proof can be carried through, apart from the application of elementary rules of vector 
analysis, with the aid of the relations: 
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(s + ^ +x2 ) r ° (xp)=45(,r)5(2/)? 

duo sin J co 6—7) V =— 7r5 r (*—-)• 



Once again, the predominant role of the lower frequencies in (22) is obvious. In fact, the 
amplitudes of the frequencies co<co^ depend numerically on an oscillating Neumann function, 
those of the higher frequencies aT>co pl on the exponentially decaying K junction. 

We mention also the following alternative representations equivalent to (22) : 



i ( * b u d 2 \ f° 
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K ° \ P \ v 2 ~? ^~'^»)' 



(24a) 



H=—cxa\ 

ire, 



{ Z z !*_"_ d«e" iu( ' f) tf„ (p^5~5 V« 2 -«?«)} • (24b) 



The square root V^ 2 - w^ is to be considered here as positive for cu 2 >co^, whereas in the 
section — o>p/<co<«^ we have to take the average of the values of the K function just above and 
just below this crosscut of the square root. The average follows from the identity: 

K U=~z)=lK (f^)+\K» (e^Vi)=~| ^o(Vi)- 

The denominator in (24a) vanishes when performing the differentiations with respect to z and t . 

5.2. The Field Expressed in Terms of Elementary Functions 

Both the electric and the magnetic field can be reduced to elementary functions by applying 
further special formulas connected with Sommerf eld's integrals. The final expressions obtained, 
viz: 



E=-(grad+5«,|) 



cos 



•fes/P-aW -" 1 } 

— — ■ . ^— i 

V( i- ?) p2+(s_<2 



H=~ curl 
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COS 






(25) 



can also be verified straightforwardly by substitution into the Maxwellian equations (18), the 
right-hand side of (18b) again being replaced by (23). The proof depends on the relation 
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{WHW-^ 2 } 



1 ^W^ 1 i\\!!!iZWIiviiri 



which results from the general relation (20) . 



=-^S(x)8(y)s(t-iy 
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According to (25) the plasma induces oscillations in the field existing otherwise in free 
space, which is in contrast with the exponential influence of the plasma in the case of a zero 
magnetic field (described in section 3.3). 

6. General Conclusions 

The described situations, referring either to a zero or to an infinite magnetic field, show 
that the field disturbances produced by a current element which is only effective at a special 
moment can, after their start at the point of observation, be decomposed into a spectrum con- 
taining exclusively low frequencies. The corresponding disturbances produced by a charge 
moving along a rectilinear path (assumed as a line of force in the case of the infinite magnetic 
field) comprise predominantly low frequencies, though theoretically all frequencies occur. 
The limitation to low frequencies is particularly pronounced at large distances from the path 
of the moving charge. These general results are expected to hold also for a finite magnetic 
field. According to preliminary calculations the spectral density then becomes infinite at 
frequencies for which the phase velocity equals that of the moving charge; therefore these 
frequencies should be observable first of all. This is in accordance with Gallet's [1959] theory 
concerning exospheric VLF emissions; in fact, the frequencies radiated from a charged particle 
moving through the exosphere are derived there from the condition expressing the equality 
of the two velocities mentioned. The general problem of a charge moving through a finite 
static magnetic field has also been approached, using a different analysis in the case of a spiral 
motion, by Pakhomov, Aleksin, and Stepanov [1962]. 
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